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Abstract 

We set up the framework for the calculation of electric dipole moments (EDMs) 
of light nuclei using the systematic expansion provided by chiral effective field theory 
(EFT). We take into account parity (P) and time-reversal (T) violation which, at 
the quark-gluon level, originates from the QCD vacuum angle and dimension-six 
operators capturing physics beyond the Standard Model. We argue that EDMs 
of light nuclei can be expressed in terms of six low-energy constants that appear 
in the P- and T-violating nuclear potential and electric current. As examples, we 
calculate the EDMs of the deuteron, the triton, and ^He in leading order in the 
EFT expansion. 



1 Introduction 



Permanent electric dipole moments (EDMs) [Tl[2] of particles, nuclei, atoms, and molecules 
provide stringent bounds on sources of parity (P) and time-reversal (T) violation beyond 
the phase of the quark-mixing matrix ^ in the Standard Model (SM). Experiments in 
preparation |1] aim to improve the current bound on the neutron EDM, \dn\ < 2.9-10"^^ e 
cm |5], by up to two orders of magnitude. At the same time, proposals [6] exist for the 
measurement of EDMs of charged particles in storage rings, in particular the proton — 
for which an indirect bound, \dp\ < 7.9 ■ 10~^^e cm, has been extracted from the ^^^Hg 
EDM [7] — and the deuteron. It may be possible, moreover, to measure in this way the 
EDMs of other light nuclei, for example ^He and ^H, as well. A nonzero EDM signal 
in this new generation of experiments would be an unambiguous sign of new physics at 
energy scales similar to those probed by the LHC, since the EDMs resulting from the 
quark-mixing matrix [8j are orders of magnitude below the current experimental limits. 

An important outstanding question is whether it will be possible to identify from the 
experimental results the microscopic parity- and time-reversal-violating (/^) source. The 
SM contains, apart from the quark-mixing phase and its lepton analog, also the fl" QCD 
6 term [9J. This interaction has dimension four and one might expect it to give the main 
contribution to hadronic f^. However, since the experimental upper limit on the neutron 
EDM constrains the vacuum angle 6 to be unnaturally small jTUl E], ^<10~^°, non- 
SM contributions from higher-dimensional/^ sources can be relevant, or even dominant. 
These higher-dimensional operators originate from an ultraviolet complete theory at a 
high energy scale beyond the electro weak scale. The first such effective operators 
one can write down have effective dimension six [121 [13 El US); the quark EDM 
(qEDM), the quark and gluon chromo-electric dipole moments (qCEDM and gCEDM, 
respectively), and two four-quark (FQ) interactions. 

It is difficult to calculate hadronic properties directly from a Lagrangian written in 
terms of quark and gluon fields. Still, long-distance effects of the strong interactions can 
be described independently of assumptions about the dynamics of QCD if one uses the 
low-energy effective field theory (EFT) known as chiral perturbation theory (xPT) [TU|[T7] 
(for reviews see Refs. [ISllin]), which translates microscopic operators into operators that 
contain nucleon, pion, and photon fields. After the translation, one is able to calculate 
hadronic properties directly from the effective Lagrangian. For the nucleon EDM (and its 
associated form factor) stemming from the dimension-four 6 term, this method was first 
used in Refs. [lOl [H], and later extended in the context of SU{2) x SU{2) [201 [2ll [221 [23] 
and SU{3) x SU{3) [21 ^ xPT. In this way it is possible to establish a hnk to PT- 
conserving (PT) interactions [211 ESI [23] . 

This approach can be generalized to include the dimension-six sources [261 [231 |27] . 
Different fl" sources at the quark-gluon level produce different hadronic interactions. The 
effective chiral Lagrangian includes not only interactions that stem from spontaneous 
chiral symmetry breaking and are therefore chiral invariant, but also interactions that 
break chiral symmetry in the same way as chiral-symmetry-breaking operators at the 
QCD level. Thus, while they all break P and T, the dimension-six operators break chiral 
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symmetry differently from each otfier and from tlie 9 term. Given enougli observables 
it sliould be possible to separate tlie various sources on tlie basis of tlie liadronic 
interactions that they generate. 

Recently it was argued [28] that a measurement of the deuteron EDM in combination 
with the neutron or proton EDM could partially separate the fundamental sources. 
A measurement of the deuteron EDM significantly larger than the nucleon EDM would 
point toward new physics in the guise of a quark chromo-EDM. The calculation was 
based on a perturbative-pion approach [291 EO] to nuclear EFT, which assumes that pion 
exchange can be treated in perturbation theory. A major advantage of this approach is 
that analytical results can be obtained. On the other hand, such a framework is applicable 
only below the momentum scale (~ 300 MeV) at which one-pion-exchange (OPE) becomes 
significant. This is the case for nuclei where the binding momentum per nucleon is small 
compared to the pion mass, but even then the size of uncertainties is set by the inverse 
of the relatively low energy scale. 

Our goal in this article is to provide a framework for the calculation of the EDMs of 
light nuclei using chiral EFT with nonperturbative OPE [ST], 15 ^ 150] . By treating OPE 
nonperturbatively, the EFT gets extended to higher momenta and thus denser nuclei, 
and convergence improves. The fact that nuclear binding momenta are small in the 
typical scale of QCD (~ 1 GeV) is sufficient for a general power counting that is able 
to estimate which hadronic interactions are dominant for each fundamental source. 
The potential has been derived previously [34j, and here we obtain the associated fj" 
currents. As explicit examples we consider the EDMs of the deuteron (^H), the triton 
(^H), and the helion (=^He) f\ 

The EDMs of the deuteron [361 EZl [381 EH HQ] and helion [HI |12] have been investi- 
gated previously within traditional meson-exchange frameworks. In the most comprehen- 
sive studies [HU] W2\ one started from "realistic" nuclear-force models and a general 
nucleon- nucleon [NN) interaction [3Sl US] . The EDMs were expressed in terms of three 
non-derivative pion-nucleon interactions, which are often assumed to be of similar 
size and dominate the EDMs, and in addition short-range interactions due to the ex- 
change of heavier mesons were included. The major advantage of a chiral EFT framework 
is that it has a direct link to QCD and exploits the chiral properties of the fundamen- 
tal sources. Moreover, the power-counting scheme allows a perturbative expansion 
such that the theoretical uncertainties can be estimated and the results can be improved 
systematically. 

When the chiral-symmetry properties of the dimension-four and dimension-six opera- 
tors are considered new insights are in fact obtained [2T1 1221 [221 [2HI At leading 
order, only two of the three pion-nucleon interactions contribute. Moreover, there are 
in general at the same order more contributions, viz. short-range contributions to the 
neutron and proton EDMs and two NN contact interactions. As we will demonstrate 
below, the EDMs of light nuclei can be expressed in terms of these six parameters, 
or low-energy constants (LECs). (Other moments, such as the deuteron magnetic 

^Similar calculations are being carried out by a group at FZ Jiilich [35j . 
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quadrupole moment, depend in addition on pion-nucleon-photon interactions [28j.) 
For three of the four sources, only a subset of these six LECs is in fact needed. Each 
LEG can in principle be calculated from the underlying source using an explicit so- 
lution of QCD at low energies, for example through lattice simulations. Compared to 
nucleons, the EDMs of light nuclei can give crucial complementary information about 
the fundamental source. However, the conventional assumption that the three 
pion-nucleon interactions can cover the whole range of nuclear EDMs is oversimplified. 

For the PT potential we use here realistic phenomenological potentials [lU HSl H6] . 
This "hybrid" approach [U] is justified whenever there is little sensitivity to the details of 
short-range physics, since such realistic potentials all include the long-range pion exchange 
that appears in chiral EFT at LO. Such an approach has been tested successfully for other 
observables [30] , such as the PT form factors of the deuteron |18] and fT NN observables 
|49j . The results in Refs. [391112] suggest that the same is true for EDMs, and we partially 
confirm this below. We use the codes from Refs. [391 112], but we recast and extend the 
results in the framework of chiral EFT with nonperturbative OPE. In particular, we apply 
power counting in order to make more model-independent statements. The cases of the 
helion and the triton are typical of a generic nucleus. However, in the deuteron, because of 
its isoscalar character, the formally LO contribution from the 9 term vanishes [3811391128] . 
a property expected [50] for nuclei with equal number of protons and neutrons, N = Z. 
We exploit the systematic character of EFT to extend the deuteron calculation for the 6 
term to the first non-vanishing order. 

Our article is organized as follows. In Section 2, we present the P- and T-conserving and 
violating interactions relevant for the calculation of light nuclear EDMs. In Section 3 we 
discuss in general the power counting of the various contributions, and present the leading 
fl" potentials and currents, while in Section 4 we specifically address nuclei with N = Z. 
Next, we evaluate the EDM of the deuteron in Section 5 and the EDMs of the helion 
and the triton in Section 6. In Section 7 we discuss our results and their implications. 
Appendices are devoted to the expression of potential and currents in coordinate space. 

2 Chiral Perturbation Theory 

xPT is the EFT of QCD for processes involving momenta Q ~ m^r <^ Mqqd, where m,r is 
the pion mass and Mqcd ~ 1 GeV is the characteristic scale of QCD. At such momenta the 
relevant degrees of freedom are nucleons, photons, and pions. The (approximate) chiral 
symmetry of QCD, SUl{2) x SUr{2) ~ 5*0(4), plays a very important role, because it 
constrains the form of the interactions involving the (pseudo) Goldstone bosons associated 
with its spontaneous breaking, the pions. In this section we review these interactions, in 
both PT and sectors of the theory. 

The xPT Lagrangian contains all interactions allowed by the symmetries of QCD. 
Each interaction is written in terms of pion (tt), nucleon (A^), and photon {A^) fields. The 
constraints of (global) chiral and (gauge) electromagnetic symmetries can be incorporated 
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through the use of covariant derivatives, 



(D^tt)^ = [d^Sab + eAf^esab) Tib (1) 



for the pion, 



N (2) 



for the nucleon, and 



= d^A, - d,A^ (3) 
for the photon. Here ~ 185 MeV is the pion decay constant, 

^-1 + 1^' (4) 

TT 

T are the Pauh matrices in isospin space, and e is the proton electric charge. Since 
~ Mqcd, nucleons are approximately non-relativistic in the processes of interest, and 
Lorentz invariance is incorporated order by order in the EFT expansion |T7]. We denote 
the (small) nucleon four- velocity by and its spin by S^; in the nucleon rest frame, 
Vfj^ = (1, 0) and = (0, c?/2) in terms of the Pauli matrices in spin space, a. 

Chiral-invariant interactions are built with pion covariant derivatives, while explicit 
chiral symmetry breaking is introduced by the average quark mass rh = (m„ + md)/^, by 
the quark mass difference md — mu = 2rhe, by electromagnetic interactions, andf and/or 
^ interactions. Each interaction with the correct symmetry transformation properties 
has a strength determined by details of the QCD dynamics. Until they are known, they 
are estimated using naive dimensional analysis (NDA) fi6[ [51], [13]. For example, the pion 
mass term originates from explicit chiral-symmetry breaking by the average quark mass, 
so = 0{rhMQc-£,). The LECs of other chiral-breaking interactions proportional to 
powers of rh and rhe can then be written in terms of m^/Mqco- Exchange of hard photons 
leads to interactions among hadrons that are proportional to the fine-structure constant 
«em = e^/47r. For simplicity, we count e ~ 1/3 as 0{1) and aem/47r as O {ml / Mqqj^) , 
since numerically aem/4vr ~ em^/(27rF^)^. It is convenient to organize the infinity of 
effective interactions in the Lagrangian using an integer "chiral index" A and the number 
/ of fermion fields [TS] : 

oo 
A=0 / 

where A = d + f/2 — 2>0, with d the number of covariant derivatives and powers of 
m^r. The index A tracks the number of powers of Mq^j-,. 

2.1 P- and T-conserving chiral Lagrangian 

The calculation of the potential and currents, which we need in order to evaluate 
nuclear EDMs, requires also PT interactions, which we summarize here. (A more complete 
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list can be found in, for example, Refs. [191 EOl |52l |32l |53].) These interactions result 
from the quark (color-gauged) kinetic and mass terms in the QCD Lagrangian. 
The terms we need in the PT Lagrangian are 



V<2,PT - 2^^" ^ 2D ^ ^ 2m 
6ml - 6m% , ^ 2\ 2 



N 



2D 



2 



(tt^ - TTg) - i^T^l - (Sm^ + 6mN^ {n x v ■ Dn) 



'^^^-D^TT ■ NrSf'N + ^ (d^tts - ^^tt ■ D^tt ] NS^N 



[iN{T X 77)3 5 ■ + H.c] 



gA^mN 



F^mN 

-TT^e'^^^'F^u {tN [1 + 2kq + (1 + 2/€i) rg] Sj:)pN + H.c.} . (6) 

The pion kinetic and mass terms have A = 0. For notational simplicity, we choose to 
absorb in the pion mass m^^ a correction oc m^. At A = 1, the leading electromag- 
netic contribution to the pion mass sphtting appears, Sm^ = 0(aem^QCD/47r), while the 
quark-mass difference contribution, 6ml = ^(^^"^^/^qcd)' is smaller by one power of 
em-rr/MQCD- The pion mass splitting, m^± — m^o = 6ml — 6ml = (35.5 MeV)^ [M], is 
dominated by the electromagnetic contribution. The nucleon kinetic terms have A = 0, 1. 
Again for simplicity, the average nucleon mass m^ absorbs a correction oc m, the sigma 
term. Entering at A = 1, 2, the nucleon mass splitting, m„ — mp = 6m + (5m jy = 1.29 
MeV |54j also receives contributions from electromagnetism and from the quark masses. 
In this case, the quark-mass contribution 5m at is expected to be the largest. By di- 
mensional analysis 6m ^ = 0{eml/MQc-D), and lattice simulations estimate it to be 
6mN = 2.26 ± 0.57 ± 0.42 ± 0.10 MeV [55], which is in agreement with an extraction 
from charge-symmetry breaking in the pn — )■ dir^ reaction [56j. The electromagnetic 
contribution is Smjy = 0(aem^QCD/4vr), which is (^(em^/MQQp) and about the 20% of 

6mN- From the Cottingham sum rule 1^T\ one finds 6mj\i = —(0.76 ± 0.30) MeV, which 
is consistent with dimensional analysis. In order to achieve the form ^ for £pT, we 
have used a field redefinition j53j to eliminate the nucleon mass difference term in favor 
of pionic mass and interaction terms. In this way, the nucleon mass to be used in nucleon 
propagators is simply m^r. The operator with LEG qa is the usual pion- nucleon axial 
coupling appearing at A = 0. We also absorb subleading corrections in it, so that the 
Goldberger-Treiman relation for the strong pion-nucleon constant, gNNn = '^gAf^N/F-,,, 
apphes without an explicit discrepancy. If for the pion-nucleon coupling constant we use 
gNNv = 13.07 [58], then gA = 1.29. Its isospin- violating counterpart is the operator with 
LEG (3i = 0{eml/M^fjY)) at A = 2. At present there are only bounds on (3i from isospin 
violation in AW scattering. The Nijmegen partial-wave analysis of NN scattering data 
gives /9i = (0 ± 9) ■ 10~^ [SH], which is comparable to estimates of /3i from ir-r] mixing. 
At A = 2 there is another isospin-violating pion-nucleon interaction generated by nu- 
cleon recoil and the nucleon mass difference. Finally, also at A = 2 there is a relativistic 
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correction |60j to the electromagnetic coupling of the nucleon involving the isoscalar and 
isovector components of the anomalous magnetic moment, respectively kq = —0.12 and 
Ki = 3.7. 

2.2 P- and T- violating chiral Lagrangian 

The lowest-dimension operator that can be added to the PT QCD Lagrangian is the 
dimension-four 6 term. With an appropriate choice of the quark fields q = {u, d)'^, the 6 
term can be expressed complex mass term [TO] . 

^rr,dim=4 = fnj qij^^q , (7) 

where = mum^/ {rriu + md) = O^ml/MQcr,) and 9 is the QCD vacuum angle, here 
already assumed to be small, 9 < 10~^°, as indicated by the experimental bound on the 
neutron EDM. 

The smallness of 9 leaves room for other sources in the strong interactions, which 
have their origin in an ultraviolet-complete theory at a high energy scale, such as, for 
example, supersymmetric extensions of the Standard Model [15]. Well below the scale 
Mr characteristic of T violation, we expect effects to be captured by the lowest- 
dimension interactions among Standard Model fields that respect SUc{S) x SUl{2) x f/y(l) 
gauge symmetry. Above Mqcd, strong interactions are described by the most general 
Lagrangian with Lorentz, color, and electromagnetic gauge invariance among the lightest 
quarks, gluons, and photons. The effectively dimension-six terms at this scale can be 
written as [HI [131 [Til US] 

CrrAim=6 = (do + d^Ts) a^'i^h F,, - ^q [do + d^r,^ a^^'^i^^'X^q G% 

^d^^,uX.fabcQa^^b^P^c^^ + ^ImSi {qqqz^'q - qrq ■ qrt^'q) 

+ ^ImS8 (gA"g qX'^i-f^q - qX" rq ■ qX^ rij^q) , (8) 

in terms of the gluon field strength Gjjj,, the Gell-Mann matrices A" in color space, and 
the associated structure constants f""^^. In Eq. ^ the first (second) term represents the 
isoscalar d^ (do) and isovector d^ (d^) components of the qEDM (qCEDM). Although these 
interactions have canonical dimension five, they originate just above the Standard Model 
scale Mw from dimension-six operators [12] involving in addition the carrier of electroweak 
symmetry breaking, the Higgs field. They are thus proportional to the vacuum expectation 
value of the Higgs field, which can be traded in for the ratio of the quark mass to Yukawa 
coupling, nig/ fq. Writing the proportionality constant as e6qfq/M^ {AirSgfg/Mj,), we have 

(ffi \ ~ I ~ ffi \ 

e5o,3^1 , rfo,3 ~ ( 47r5o,3^ 1 , (9) 
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in terms of the average light-quark mass rh and the dimensionless factors 60,3 and 5o,3 that 
represent typical values of 6q and 6q. The third term in Eq. ^ [I3] is the gCEDM, with 
coefficient 

in terms of a dimensionless parameter w. The fourth and fifth operators [TH [15] are 
FQ operators, with coefficients 

in terms of further dimensionless parameters ai^g. The sizes of (5o,3, <5o,3; and di s 
depend on the exact mechanisms of electroweak and PT breaking and on the running to 
low energies where nonperturbative QCD sets in. The minimal assumption is that they 
are C(l), 0{gs/4:7i), C(((7s/47r)^), and C(l), respectively, with gs the strong coupling 
constant. However, they could be significantly smaller, when parameters encoding 
beyond the Standard Model are small, or significantly larger, since fg is unnaturally small; 
for discussion and examples, see for instance Refs. [21 [15]. 

The operators in Eqs. ([T]) and ([S]) have different transformation properties under chiral 
symmetry, which has profound implications for the form and relative importance of the 
fl" pion-nucleon and NN couplings in the effective Lagrangian. The 9 term in Eq. (j7| 
transforms as the fourth component of an S'0(4) vector P = {qTq,qi'y5q), the third 
component of which is responsible for quark-mass isospin violation [32] • from the 6 
term and isospin violation from the quark mass difference are therefore intrinsically linked; 
this link appears in certain relations [TOl Ull [22] between the coefficients of and isospin- 
breaking operators in xPT through a coefficient p = (1 — e^)d/2e. The dimension-six 
operators in Eq. ^ have different transformation properties still [2Zl [221 1231 [2H1 El] ■ The 
isoscalar and isovector qEDM and qCEDM transform as the fourth and third components 
of two other 5*0 (4) vectors. There is no useful link to PT observables, and the third 
component of the qCEDM vector tends to generate hadronic interactions, which for 9 
require tensor products and are of higher order. For qEDM, purely hadronic interactions 
arise from integrating out at least one hard photon, which leads to further breaking of 
chiral symmetry in the form of tensor products of the vectors with an antisymmetric chiral 
tensor [32] . The contributions of the qEDM to hadronic couplings, like pion-nucleon or 
NN couplings, are suppressed by aem/47r. In contrast, the gCEDM and the two fj" FQ 
operators are singlets of the chiral group. Because they are chiral invariant, and contain 
no photons, the gCEDM and the two FQ operators lead to exactly the same effective 
interactions, although, of course, with different strengths. For simplicity of notation, 
in the following we treat gCEDM and FQ operators together; we refer to them as 
chiral-invariant (%!) sources and use w to denote both w and aifi. 

{w,ai,as}^w. (12) 
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We now present a subset of the complete chiral Lagrangian originating from the 
fundamental sources above. We only give the operators that play a role in the LO calcu- 
lation of light-nuclei EDMs, the more general Lagrangian being found in Refs. [221 Ell- 
in general, a LO calculation of the EDM of a light nucleus requires six interactions: 

C„ = -2N{do + d,T,)S^Nv''F^,-^N{goT-7r + g,rr,)N 

+CiNNdf,{NS^'N) + C2NTN ■d^{NS^TN) + ... , (13) 

which represent short-range isoscalar (do) and isovector (Ji) contributions to the nucleon 
EDM, isoscalar (go) and isovector (^i) non-derivative pion-nucleon couplings, and two 
short-range AW interactions (Ci, C2). Here we relegate to the ". . ." terms related to 
the above by chiral symmetry. The explicit forms of these terms depend on the fl" source 
but, because they involve more pion fields, they do not appear in the LO EDMs we are 



interested in. Note that Eq. (13) is the form of Cfj- after a field redefinition is performed 
to eliminate pion tadpoles and guarantee vacuum alignment; the parameters thus absorb 
contributions generated by this field redefinition. 

Which of these six interactions is relevant depends on the system we are studying and 
on the fundamental source. As will be seen, the spin and isospin of the deuteron cause 
the deuteron EDM to be sensitive to only three of the above operators. In more general 
cases, the EDMs of light nuclei are sensitive to all six interactions. The EDMs of heavy 
nuclei could involve more operators than the set above. Generically one might expect a 
dominance by effects from (z) a single nucleon, since multi-nucleon contributions tend to 
be suppressed at low energies by phase space; and {ii) pions, thanks to their small mass 
and related long range. However, significant deviation from this expectation comes from 
the relative sizes of the various LECs, which depends on the fl" source. NDA leads to 
the following estimates for the dimension-four and -six fl" sources: 

• For the 6 term, four operators play a role at LO, the other two appearing only 
at sub leading orders. In order to generate gi, which is relevant for the deuteron 
EDM, the 6 term requires an insertion of the quark mass difference, which causes 
a suppression of gi relative to go by a factor eml / Mq^j) [22]. (At the same order, 
there exists also a two-derivative pion-nucleon coupling, but for our purpose here it 
can be absorbed by a small change in g^ ^].) The LECs scale as 



^-K^"^)' '"^''ri'- 

For the qCEDM, the same four operators are needed. In this case, there is no a 
priori relative suppression of gi and the LECs scale as 

do, - ^1(^0 + ^3)^^). (15) 
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(Here the "+" signs are not to be taken literally; they are only meant to signify two 
independent contributions to a LEG.) 



For the qEDM, only the short-range EDM contributions are important, and they 
scale as 

For the xl (gCEDM and FQ) sources, the non- derivative pion-nucleon interac- 
tions, which break chiral symmetry, are suppressed by a factor m^/Mq^j-, compared 
to short-range nucleon EDM contributions and NN interactions, which conserve 
chiral symmetry. (Again, a two-derivative pion-nucleon interaction exists at the 
same order but can be absorbed in qq [34J.) All six operators thus become relevant, 
and the LECs scale as 




2.3 EDM of the nucleon 

Using these interactions, the nucleon EDM has been calculated in xPT up to NLO for all 
sources of dimension up to six [251 ES, [23] . 

In the power counting of xPT [16], one considers typical momenta Q ~ ~ ^ 
^QCD ^ rriN 271 Fj, and assigns 

• a factor Q^/(47r)^ for each loop integral; 

• a factor 1/Q for each nucleon propagator; 

• a factor l/Q"^ for each pion propagator; 

• the NDA estimate for the LECs corresponding to the interactions in the diagram. 

This produces for any observable an expansion in the small ratio Q/Mqcd- 

For the nucleon EDM, in all cases there are short-range contributions from do^i at LO. 
For qEDM and xl sources, the relative suppression of pion-nucleon couplings means that 
loops come at higher orders and only (Jo,i appear up to NNLO [2^. In contrast, for 9 and 
qCEDM, one-loop diagrams contribute at LO and NLO. Using dimensional regularization 
in d dimensions at a renormalization scale /i, and introducing 

rT _ egAgo / 2 47r/i2\ 
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where 7^ ~ 0.577 is the Euler-Mascheroni constant, the isoscalar and isovector EDMs 
can be expressed at NLO respectively as [25], EH [23] 



Sm^r Sm 



Am 



N 



N 



+ 



and 



di = di + 6di + 



egAgo 



In 



N 



+ 



57r m,r 
4 miv 



mt 



+ 



(19) 



(20) 



where the gi terms apphes to qCEDM only. The dependence on the arbitrary scale /i in 
5di is compensated by di. In fact, the loop contributions cannot be separated from the 
short-range pieces in a model-independent way. After absorbing all these terms in rfo,i! 
which we do for the rest of the paper, we can write for all sources 



for the neutron and 



dn — do — di 



dp = dQ + di 



(21) 
(22) 



for the proton. 

However, one expects no cancellation between short-range contributions, which are 
analytic in m^, and the "chiral-log" and other finite terms, which are not. Thus the 
non-analytic terms serve as lower-bound estimates for the size of dp^n- We then expect, 
for e dl] and qCEDM 



do > 0.01 
dl ~ 0.1 



^ + 0.3^ 
F F 

1^ + 0.03|^ 
F F 



efm 



efm 



(23) 
(24) 



As stressed in Ref . [26] , measurements of both dn and dp alone can tell us little about 
the underlying source of f^. More can be learned from measuring the EDMs of light 
nuclei, the calculation of which we now turn to. 



3 Ingredients: the generic case 

The EDM of a nucleus with A> 2 nucleons can be separated into two contributions. The 
first contribution comes from an insertion of the electromagnetic current J^. The 
second stems from the PT charge density Jprp upon perturbing the wavefunction of the 
nucleus with the potential Vf^, such that the wavefunction obtains a component. 
To first order in the fl" sources, the EDM is thus a sum of two reduced matrix elements 



dA 



D 



PT 



(25) 
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The nuclear ground state |^I/a) and its parity admixture {"^a) are the solutions of homo- 
geneous and inhomogeneous Schrodinger equations, 

(E-ifpT)I^A) = 0, (26) 
{E-Hpt)\^a) = V„\^a) , (27) 

respectively, where HpT is the PT Hamiltonian. The potential Vf^ is shown in 
coordinate space in Appendix A The EDM operators DpT and D^r^ are obtained from 
the corresponding charge densities Jpj. and J^, respectively, as discussed in Appendix 
[B| The factor of 2 in front of the second matrix element corresponds to the number of 
time-ordered diagrams, and the phases of wavefunctions are chosen so that these matrix 
elements are purely real. 

In this section we identify the ingredients needed for the LO calculation of (Ia, assuming 
no particular cancellations or suppressions due to spin/isospin factors. 



3.1 Power counting 

Both the potential and the current can be obtained from the Lagrangian of the 
previous section. The potential Vf^ for the various sources has been derived in Ref. 
|34j . To the order we are concerned with here, the potential can be taken as two-body. 
The and PT currents can also be divided into one-body and more-body currents. 
As we will see, the latter are dominated by two-body effects as well. There are thus 
four classes of contributions to a nuclear EDM, schematically drawn in Fig. [T} In order 
to determine which diagram(s) give(s) the most important contribution(s) we need to 
estimate their sizes by applying power counting. 

We need to count powers of the generic momentum Q in the process, in order to get an 
expansion in Q/Mqcd- Here Q is given by the nuclear binding momentum, which for a 
typical nucleus can be taken as Q ~ Tn.,^ ~ F^r, as standard in xPT. However, as pointed 
out by Weinberg pi], the power counting of xPT needs to be adapted to the existence 
for A > 2 of intermediate states consisting purely of propagating nucleons. A generic 
diagram can be split into "reducible" parts, that contain such states, and "irreducible" 
subdiagrams, which do not. Within an irreducible subloop, the contour integration over 
the 0th component of the loop momentum can always be performed in such a way as to 
avoid the nucleon pole. In these diagrams the nucleon energy is of order Q, as assumed 
in xPT power counting. On the other hand, in diagrams where the intermediate state 
consists purely of propagating nucleons, i.e. reducible diagrams, one cannot avoid the 
poles of nucleon propagators, thus picking up energies ~ Q'^/tiin [31] rather than ~ Q. 
Moreover, such loops also obtain an additional enhancement of An. The contribution of 
such a reducible diagram can be counted by applying the modified rules [30j : 

• a factor /{iTrniN) for each loop integral; 

• a factor ttin/Q^ for each nucleon propagator; 

• a factor l/Q"^ for each pion propagator; 
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a bed 




Figure 1: The four general classes of diagrams contributing to a nuclear EDM described 
in the text. Solid and wavy lines represent nucleons and photons. The three (two) dots 
stand for A — 3 {A — 4) nucleon propagators. The large triangle denotes the nuclear 
wavefunction; the oval, iterations of the PT potential; the dot with an attached photon, 
the PT one-body current; the oval with an attached photon, the PT two-body current; 
the black square, the potential; and the black square with an attached photon, the 
fl" current. 



• the NDA estimate for the LECs corresponding to the interactions in the diagram. 

As an example, consider an insertion of a LO, PT pion exchange in a diagram. It gives 
rise to one additional loop ~ /{inniN), two nucleon propagators ~ mj^/Q^, a pion 
propagator ~ and two insertions of the strong pion-nucleon vertex ~ {Q/F-^y. 

Combining these factors, the extra one-pion exchange amounts to Q/M]\f^, where M^jy = 
AirF^/mj^r ~ F^^. A similar power counting holds for short-range PT interactions, although 
the situation for them is more complicated [33] • For very light nuclei, Q < M^^ and pion 
exchange can be treated perturbatively [291 EH] • The deuteron EDM has in fact already 
been considered in this light [26]. For less dilute nuclei, however, one expects Q ~ M^n 
and pion exchange needs to be summed to all orders [301 ES] • The counting rules above 
are a generalization for A > 2 of the rules given in Ref. [30]- Note that they provide 
a natural explanation for the Q/Mqqd supression associated with an additional nucleon 
observed in pion-nucleus scattering [6ll [35| . 

We can now estimate the size of each of the classes of diagrams in Fig. [Tj For each 
class we take the PT and LO interactions in Eqs. ^ and (13), respectively. The 



iteration of the LO PT potential costs no factors, and is necessary among nucleons in 
reducible intermediate states, as indicated in diagrams (c) and (d) of Fig. [l] Such iteration 
among nucleons before and after all fl" and electromagnetic insertions builds up the PT 
wavefunction, represented in Fig. [Tj as well, which introduces an overall normalization of 
the diagrams. This normalization can be read off from the diagram analogous to (a), where 
the one-body current is given instead by the electromagnetic charge. In the following we 
account for this normalization by omitting the A — 1 loops and A + 1 nucleon propagators 
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that are common to all diagrams. Thus, diagram (a) is simply 



Da = (dp^nQ) ■ 



{2i 



In contrast, diagram (b) has one additional irreducible loop ~ / (47rm7v), one additional 
nucleon propagator ~ vn^jQ'^^ and the leading two-body current. For both qEDM 
and y\ sources the latter brings a suppression of a factor Q^/Mqqq, whereas for the other 
sources the contribution is comparable to the one-body term. One can continue in this 
fashion to find that for diagram (c). 



9o,i ^ 



(29) 



while for diagram (d) there is always a further suppression by a factor Q^/Mq^j-,. Anal- 
ogously, more-body potentials and currents bring further suppression. 



Plugging in the scaling of the LECs for the different sources, Eqs. (14), (15), (16), and 



(17), we can draw the following general expectations for the EDMs of light nuclei: 



• For the 6 term, the nuclear EDM is dominated by diagram (c): the nuclear wave- 
function acquires a admixture after a one-pion exchange involving the isoscalar 
go vertex; the admixed wavefunction then couples to the proton charge. 

• For the qCEDM, the nuclear EDM is dominated by the same effect as the 6 term. 
However, for the qCEDM the fl" pion- nucleon vertex can be either qq or gi. 

• For the qEDM, the nuclear EDM is dominated by the sum of the EDMs of the 
constituent nucleons, diagram (a). 

• For xl sources, the nuclear EDM is more complicated than for the other sources. 
Due to the chiral suppression of the pion-nucleon interactions, diagrams (a) and (c) 
are equally important, and in the latter the short-range fl" interactions (7i 2 need 
to be included besides the one-pion exchange from both go and gi couplings. 

3.2 P- and T-odd potential 

For all sources considered, except qEDM, an insertion of the LO two-nucleon potential 
appears in the EDM at LO. The general NN potential was derived in Ref. [3lj and 
we summarize the relevant parts here. In momentum space the potential is given by 



k 



+1 



F2 

9Agi 



2F3 



k'^ + ml 



-(0 



To — T, 



(i) 



k"^ + ml 
(30) 
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where /2 (r^") /2) is the spin (isospin) vector of nucleon n, and k = pi — is the 
momentum transferred from nucleon i. In this expression, at LO go originates from 6- 
term, qCEDM, and xl sources; gi from qCEDM and sources; and Ci from sources 
only. The pion-exchange parts are well known (for example, Refs. [501 EHl ES]), while 
the contact interactions incorporate all other effects of short-range, such as single 
exchanges of the mesons u and r] (Ci) and p (6*2) [Mj- 

3.3 Currents 

As we argued above, only one-body currents are necessary at LO. For the 6 term, qCEDM, 
and xl sources we need the PT current coming from the proton charge in Eq. (|6]), 

4t = i (1 + ' ^^^^ 

where T^^y2 is the isospin of the nucleon that couples to the one-body current. 

For the qEDM and xl sources we need as well the current originating from the 
nucleon EDMs, 

= -^ {do + Jir«) (T« ■ q , (32) 

where a^^y2 is the spin of the nucleon that interacts with the photon and q is the outgoing 
photon momentum. 



4 Ingredients: nuclei with N = Z 

Although the power counting discussed above holds for general light nuclei, it is possible 
that a diagram, which is expected to be LO, does not contribute to the EDMs of certain 
systems. For nuclei of equal neutron and proton number, N = Z, i.e. the third component 
of isospin = 0, an insertion of the isoscalar fl" potential in combination with the LO 



one-body PT current, i.e. Eq. (31), does not contribute to the EDM ^50j. To see this, 
consider the EDM operator resulting from the LO one-body PT current, which takes the 
simple expression 

m=it{^+-^")('=li:^'''^^' (33) 

i=l i=l 

in intrinsic coordinates with XliLi 'Ci = 0. Since this operator is isovector, i.e. A/ = 
1, and conserves I3, i.e. A/3 = 0, it can only yield a non-vanishing moment when 
the nuclear state of a (/ , Is = 0) nucleus acquires some parity admixture with isospin 
(J = J ± 1 , J3 = 0). Therefore, one needs isovector components in Vfrp to induce such 
admixture. The above argument holds in the non-relativistic limit. 

This observation is of no concern for sources where there are other contributions at 
the same order as those contributions that vanish. The nuclear EDM is then simply 
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dominated by the non- vanishing LO terms. For the 9 term, however, the LO contribution 
consists only of an insertion of the isoscalar potential, such that, for N = Z nuclei, 
we need to go further down in power counting to find the dominant EDM contributions. 



4.1 Power counting 

Because the formally leading diagram (c) of Fig. [l] vanishes for N = Z in the ^-term 
case when both the PT one-body current and the two-body potential are used, let 
us first consider corrections in this diagram. It turns out that NLO corrections to both 
the fl" potential [Mj and PT one-body current vanish, and the first corrections we need 
to account for are at NNLO. By looking at the scaling of the LECs for the 6 term in 
Eq. ( 14 ) and the power counting for the classes of diagrams in Fig. [l| we then conclude 



that the first non-vanishing contributions can come from all classes of diagrams: the LO 
nucleon EDMs in diagram (a), the LO fl" two-body currents in diagram (b), the NNLO 
two-body potential or the NNLO PT one-body current in diagram (c), and the LO 
PT two-body currents with the LO two-body potential in diagram (d). 

For the other sources only parts of the LO contributions given in the previous section 
remain. For qCEDM and sources we need the potential from gi OPE. For qEDM 
and xl sources we also need the isoscalar short-range contribution to the nucleon EDM. 



4.2 P- and T-odd potential 



For qCEDM and xl sources we can use the same potential as in the generic case, but 
the go and Ci^2 terms will not contribute. We do not require a /^-potential for qEDM. 
For the 6 term we need the NNLO fl" potential calculated in Ref. [34j- At this order 
further isoscalar terms appear, which also will not contribute. Thus we need here only 
the following terms: 



V„{k,K,P) 



2F^ 

+ (^9Agi 

Ma 



9a9i - 



4' + r«) 



0" 



a 
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To — r. 
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3F2 



5rni — 5ml 



— dmi 



k 



^2 _|_ 



N 



X [a 



a 



) 



k^ + ml 
k , go9ASmM 



X 



(<?W + <jO)).ie+(a«-<T(^)) 



X T 



P ^{P-k)k 
2 k"^ + ml 



1 



k'^ + m: 



where P = pi + pj is the center-of-mass (CM) momentum of the nucleon pair and K = 
{pi +p^ — Pj —pj)/A. The first two terms originate in one-pion exchange with gi instead 
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of go or with /3i instead of gA- The next term arises from isospin breaking in the pion and 
nucleon masses, and it is very small [3l]. The last term is due to isospin breaking in the 
pion-nucleon vertex. The potential also includes corrections [31], which we do not 

include here for the reasons given below. 



4.3 Currents 

For the same reasons that require the NNLO potential we also need the NNLO PT one- 
body electric current, to be used with the ^-term LO potential. Again we do not bother 
with terms that give a vanishing contribution for = Z nuclei. The only remaining 
correction from Eq. ^ is given by 



7° 

Jprp 



%e 



{Pi+P'iY 



1 + 2^0 + (1 + 2ki) t^^ 



(35) 



which agrees with Ref. [HD]- Here Pi (p-) is the momentum of the nucleon that couples 
to the photon before (after) interaction. 

We also need two-body currents, both PT and f^. We use incoming momenta pi = 
P/2 + p and pj = P/2 — p and outgoing momenta = P '/2 + p' and pj = P' /2 — p' . 

The photon momentum q = P — Pis outgoing. For convenience we introduce k = p — p' 
as before, K = [p + p')/2, and Pt = {P + P )/2. In the evaluation of the currents at 
the order we are interested we can use the nucleon on-shell relation p° = p^/2mN, or 
alternatively = {Pt ■ k — q- K)/2mN. 

The relevant diagrams for the LO two-body PT electric current, used again in combi- 
nation with the LO two-body potential, are shown in Fig. |2} All interactions come 
from the PT Lagrangian, Eq. (|6|). In momentum space the current reads 



'JpT,a 



2ieg 



[cf«-(fc + g72)][a(^-)-(fc-g/2)] 
[{k + q/2Y + ml] [ik - q/2y + ml] 



J 





PT,b 



iegl 



2F^mN 



X T 



(i) 



X 



■ {Pt + 2K)] [<t(^-) . (fc - g/2)] ^ {Pt - 2K)] [c?W ■ {k + q/2)] 



{k - q/2y + ml 



2eg', 



PT,c 



F3 



/ q / q 



{k + g/2)2 + ml 
W-(fc + g72)][a(^-)-(fc-g72)] 
[{k + ql2Y + m2] \{k - q/2Y + ml] 



(36) 



We also need to include the LO two-nucleon electric current. The diagrams con- 
tributing to this current are shown in Fig. [3j Here PT interactions come from the PT 
Lagrangian, Eq. ([6]), and the interaction is the g^ vertex in the Lagrangian, Eq. 
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a 



b 



c 



Figure 2: Diagrams contributing to the PT two-nucleon electric current. Solid, dashed, 
and wavy lines represent nucleons, pions, and photons. A diamond marks an isospin- 
breaking PT interaction and the other vertices isospin-conserving PT interactions: lead- 
ing (filled circles) and subleading (circled circles). Only one topology per diagram is 
shown. 



Figure 3: Diagrams contributing to the fl" two-nucleon electric current. A square marks 
a interaction; other notation as in Fig. |2| Only one topology per diagram is shown. 



(13). The current is given by 



r 



TT 

__e5U|o_ 
2iegAgo 



(gW + a(j')) ■ q/2 + (a(^) - a(j)) ■ k 
' " [{k + q/2Y + ml] [{k - q/2Y + ml] 
■ {Pt + 2K) (t(^') ■ {Pt - 2K) 



5mN ( T 



{k - g/2)2 + ml {k + q/2f + ml 



(37) 



^3 ^3 



[(A; + g72)2 + ml] [{k - q/2f + ml 



5 EDM of the Deuteron 

We are now in position to calculate the EDM of the deuteron, which provides the simplest 
example of an = Z nucleus. The ground state of the deuteron is mainly a ^Si state. 
The deuteron obtains a ^Pi component after a Qq pion exchange or an insertion of Ci^2- 
Since the LO PT one-nucleon current is spin independent, it cannot bring the deuteron 
wavefunction from ^Pi to ^5*1, and therefore these contributions vanish for the deuteron, 
as anticipated on more general grounds in the previous section. 

The deuteron EDM has been studied before in the meson-exchange picture [SUj [371 
ESI EHl SO], with various degrees of sophistication in the treatments of the P- and T- 
conserving interaction HpT- Using modern high-quality phenomenological potentials 
[m Il5] , Ref . [39^ found that the model dependence of HpT is rather small for a deuteron 
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EDM generated by the OPE sector of the interaction. The detailed study in Ref. ^0] 
confirmed this point. Since our new EFT scheme shows that the leading-order contribu- 
tion from various sources to the deuteron EDM also comes from the long-range terms 
in Vfp, we take advantage of the existing calculation scheme of Ref. ^39j to obtain wave 
functions I^E^an) and |\E'2h). The calculation is performed in coordinate space using ihef^ 
potentials and currents from Appendices |X] and |B} respectively. Of course, a fully con- 
sistent treatment would involve using the PT interaction Hpx derived from the complete 
chiral Lagrangian, instead of a phenomenological potential. At present, unfortunately, 
such a consistent potential does not exist beyond LO [33]. It would include relativistic 
corrections as well, which are absent in the phenomenological potentials we use. For this 
reason, we neglect relativistic corrections in the potential and currents as well. We 
expect that the results from a fully consistent calculation will not deviate significantly 
from the results we obtain here. The numbers below correspond to the Argonne fig po- 
tential [15], but results for the Reid93 and Nijmegen II potentials agree within 5%. 
This is less than the error of order mTr/MqcD ~ 20% intrinsic to xPT in lowest order. 

The simplest contribution to the deuteron EDM comes from the constituent EDMs. 
The LO J^, given in Eq. (32), yields a one-body EDM operator 



i=l 



{31 



For the deuteron, an isoscalar (/ = 0) and spin-triplet {S = 1) state, one simply gets 



^2 



^(1) 



"$2 



(39) 



In order for -Dp-),, a purely isovector operator as discussed earlier, to yield a non-zero 



contribution in the deuteron, it is obvious that the parity admixture |\E'2jj) has to be a 
^Pi state. Among the various terms in the LO potential, Eq. (30), only the one with 
the isospin-spin operator {t^^ — T^'^){a'^^'> + a'^'^^) can contribute. The result is 
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PT 



^2h(=^Pi; 



-0.19 — efm 



(40) 



However, when it comes to the 9 term, because the LO contribution vanishes as argued in 
the previous section, the leading contribution is in fact NNLO. Among the higher-order 



interactions identified in Section 4.2, the terms with coupling constants {qaQi + 5'o/3i/2) 
and gAdoSmN in Eq. (34) can contribute, by isospin and spin selection rules. Except for 



the coupling constants, the operator structures of the former are the same as the one in 



Eq. (30), so the matrix element can simply be obtained by replacing 



9i 9i 



A 

2gA 



90 



(41) 
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in Eq. (40). Combining this with the contribution from the isospin-breaking pion-nucleon 
vertex, we find the matrix element 



\l/2 



D 



(1) 

PT 



0.19 



9i_ _^ Pi 9o 
Ftt '^9 a Ft, 



+ 5.8 ■ 10 



-4 9o 



efm (42) 



for the 9 term. 

For the 9 term, there are in addition NNLO currents to be taken into account. For the 
currents, as the corresponding EDM operators are sandwiched between two isoscalar 
states, they must be isoscalar to contribute. Among the NNLO currents identified in 
Section 4.3, only the third current in Eq. rt38), J^^' iiisets the requirement and leads to 
a two-body EDM operator (see Appendix jBT 
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(43) 

in terms of the positions Xi and x^ of the two nucleons and the derivatives Vi and V2 
with respect to them. This results in the matrix element 
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7! 



^2H 
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^2H 



1.1 ■ 10-^ — efm 



(44) 



for the deuteron EDM. The contributions of two-body PT currents to the EDM have 
again to be coupled with the parity admixture generated by the LO Vfp, which is purely 
isoscalar when 9 is the source. The only PT current with an isoscalar component, 
among those identified in Section 4.3, is the third current in Eq. (36), JpTc- gives a 
two-body EDM operator 



^(2) 



^(1) ^(2) 



(aW-V0(c?(')-V2), (fi+X2) 



g-m^|xi-X2| 



5) 

Since the isoscalar parity admixture |\I/2h) can only be a ^Pi state, this current gives a 
matrix element 
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^2H 
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^2h(^Pi; 



-3.3- 10-^1^ efm 



(46) 



In total the deuteron EDM can be written as a function of three LECs, 



d2ii = dp + dn + 



-0.19 ^ + (0.2 - 0.7 ■ 10^ /3i) ■ IQ- 



9o 



efm , 



(47) 



where dp^n should be included for 9, qEDM, and xl; 9i for 9, qCEDM, and xl; and go for 
9 only. 

This result can be compared, for each of the sources, with the calculation where OPE 
is treated perturbatively |28]. For both qCEDM and qEDM the nonperturbative pion 
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approach adopted here agrees very well with the perturbative calculation. In the case of 
the qCEDM, it was also found that the deuteron EDM is dominated by gi pion exchange 
and given by [28j 



d2H(qCEDM)| 



pert 



671 F^m^ (1 + 



-0.23 1^ efm 



(4J 



where 7 ~ 45 MeV is the binding momentum of the deuteron. This result agrees exactly 
with a zero-range model [38] and is 22% larger than the result from the qCEDM calculation 
with nonperturbative OPE [39] reproduced above, 



d2H(qCEDM) 



-0.19 1^ efm 



(49) 



Since the estimated error in the perturbative calculation is of order Q/Mtvat ~ 30%, the 
calculations agree within their uncertainty. By power counting the contribution from do 



is expected to be suppressed by m'^/MQQj) compared to Eq. (49). From Eq. (15) we infer 
doF^^/gi = C(eF^/MQQ]-,) ~ 0.03 efm, implying that, in the case of qCEDM, the nucleon 
EDMs contribute at the 30% level to the deuteron EDM. This suppression is less than 
formally expected. If we assume the isoscalar nucleon EDM is saturated by its long-range 
part, Eq. (23), the contribution is at the 10% level. In any case, the correction by the 



isoscalar nucleon EDM is of the order of the intrinsic xPT uncertainty m^/Mqco, such 
that for the qCEDM the deuteron EDM at LO is given by Eq. ( [49] ). 
Likewise, for qEDM the conclusions of Ref. 



do not change once we treat OPE 
nonperturbatively. The deuteron EDM is in this case simply the sum of the neutron and 
proton EDM, 

rf2H(qEDM) = 2do . (50) 

The comparison is more subtle for 6 and xl sources. For both of these sources, 
the deuteron EDM is expected in the perturbative-pion approach to be dominated by 
the isoscalar nucleon EDM, since pion exchange is further suppressed in the Q/M^n 
expansion. In the nonperturbative power counting/^ pion exchange is a dominant effect 
as well. In order to compare the two effects — nucleon EDMs and pion exchange — in the 
nonperturbative calculation we can look at the estimated scaling of the LECs. For xl 
sources, 



d2H(xI) = 2rfo -0.19 1^ efm . 

-TV 

OieFjeml) 



(51) 



From Eq. (17) we infer that F-j^do/gi = O^cFt^/ em^) ~ 5 efm. Thus, although formally 
gi exchange is LO, because of a combination of e suppression and the relatively small 
factor of 0.19 in Eq. (47), it actually is expected to contribute only at the ~ 5% level to 



the deuteron EDM. For 6 there are additional contributions from g^, 



d2^{e) = 2do + 



0.19 1^ + (0.2 -0.7- lO^/Ji) ■ 10 
Fn 



-3^ 



efm . 



(52) 



The contributions from the and PT two-body currents, Eqs. (43 ) and (45 ) respectively, 



are of similar size. The NN data constraint 1591 on (3i shows that the contribution from 
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the potential is no larger, and the full go term is < 0.9 ■ lQ~^{go/ F-j^) efm. From Eq. 
(14) we expect that gi/go = 0{eml/MqQY)) ~ 10^^, so the gi/F.„ contribution should be 
comparable to these small (7o/-^7r contributions. In contrast, we expect a larger weight 
from the pion cloud around each nucleon, which for (Iq enters at NLO and gives Eq. 
(23). Thus again, although pion-exchange contributions in the potential and currents are 
formally LO, e suppression and relatively small numerical factors in the deuteron make 
them likely no more than ~ 10% of the nucleon EDM contribution. 

The fact that pion-exchange contributions are expected to be smaller in the deuteron 
than assumed in xPT power counting confirms that the power counting of Ref. [28] , 
where pion exchange comes in at NLO, works better for a loosely bound nucleus. The 
XPT power counting should become more accurate as we consider heavier, denser nuclei, 
the simplest of which we tackle next. 



6 EDM of the Helion and the Triton 

In this section we investigate the EDMs of ^He and ^H. No particular cancellations are 
expected, so the framework of Section [3] applies. 

The EDM of ^He was studied in Ref. [32] , where two mechanisms were considered: 
nucleon EDMs and a two-nucleon potential containing the most general non-derivative, 
single vr-, p-, and w-meson exchanges. The nuclear wavefunction was calculated with 
the no-core shell model (NCSM) [62], where a PT nuclear potential is solved within a 
model space made from appropriately symmetrized combinations |63j of Nmax harmonic- 
oscillator wavefunctions of frequency VL. In Ref. [32] both Argonne f 18 [IS] and EFT- 
inspired [IB] potentials, including the Coulomb interaction, were used. At large enough 
Nmax results become independent of Vt. 

Here we adapt this calculation to the ingredients from chiral EFT, and calculate 
the EDM of for the first time. As argued in Section [sj power counting for generic 
light nuclei tells us that for all sources of dimension up to six, the EDM is indeed 
expected to come mostly from the nucleon EDM and from the two-nucleon potential. 



as assumed in Ref. [32]. The only difference is that the EFT potential (30) contains, 
in addition to OPE, also two LECs {Ci and C2) representing shorter-range interactions. 
This potential in coordinate space is given in Appendix |Xj The OPE terms were included 
in Ref. 02], while C*i and C2 can be thought of as originating from, respectively, u and 
p exchanges, also considered there. The relation can be made quite explicit if we choose 
to regularize the delta functions with Yukawa functions, following a strategy successfully 
employed before to study the effects of the EFT f'T potential [31] : 

'"'^'-e-'"^'- -> C'i5(3)(r-') , (53) 
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Attt 



^.m,r ^ C26^^\r) , (54) 



as mi, 2 — ^ 00. When mi = rriui (m2 = m,p) and Ci (C2) is an appropriate combination of 
ui [p) couplings [33] , the expressions on the left-hand side coincide with those in Ref. [32] • 
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Figure 4: Magnitude of the tri-nucleon EDMs in units of F^CiC fm, as function of the 
regulator mass in GeV: i = 1 (left panel) and i = 2 (right panel). The solid (dashed) 
curve is for helion (triton). 

Here we recalculate these contributions for values of mi_2 up to 2.5 GeV. For uniformity 
with Section|5]we again display numbers obtained with the Argonne f 18 potential. In Ref. 
|42j it was found that for helion the contributions from nucleon EDMs ((io,i) and from 
pion exchange (^o,i) change with PT potential by no more than ~ 25%. We have verified 
that the same is true for triton. Unfortunately the situation is different for the short-range 
two-body contributions ((5i,2), which are much more sensitive to the PT potential, as we 
discuss shortly. 

The nucleon EDM contributions are found to be 

^3He) = 0.88dn-0M7dp , (55) 
^3H^ = -0.050 d„ + 0.90 dp . (56) 

As expected, the helion (triton) EDM is mostly sensitive to the neutron (proton) EDM 

For the contribution from the potential, our results for triton are very similar in 
magnitude to those for helion, in the case of OPE already obtained in Ref. |12]. The 
contribution of Ci^2 as a function of mi 2 is given in Fig. |4]for Argonne vl8. For each 
regulator mass, we perform calculations at four values of i7 = 20,30,40,50 MeV, up 
to Nmax = 50. We observe convergence and estimate a 10% error from the spread of 
results with Q. (See Fig. 1 of Ref. [12] for a generic convergence pattern.) As it can 
be seen from Fig. [4| the results become approximately mi^2 independent at large masses, 
implying that Ci^2 approach constants in this limit. Results are very different for the 
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^3H 
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EFT-inspired potential. Within the region of masses studied, we found an approximately 
linear dependence on the regulator mass, always larger in magnitude than for Argonne 
f 18. While for mi = rrii^ and 1712 = rrip the contributions to the tri-nucleon EDMs differ 
by a factor ~ 2 [32], the difference first increases and then decreases as 2 increases, 
but it is still a factor of ~ 5 at 2.5 GeV. The linear regulator dependence could indicate a 
different running of C*! 2, or simply a very slow convergence. However, calculations with 
this potential are computationally more intensive and we have been limited to Nmax = 40, 
which increases the error. In any case, there is clearly a much stronger dependence of 
these short-range contributions on the potential, and more solid numbers have to await a 
fully consistent calculation. We quote here the Argonne vl8 numbers at 2.5 GeV, but we 
emphasize that they represent only an order of magnitude estimate. We obtain 



He 



^3H 



D 
D 



(1) 



PT 

(1) 
PT 



^3He 
^3H 



-0.15|^ - 0.28^ - OMF^Ci + 0.02F^C2 ) efm(,57) 

TV -^77 

0.15^ - 0.28^ + O.OlF^Ci - 0.02^^3^2 ) efm . (58) 



In total, then, as anticipated in Sections 2.2 andkll the EDMs of helion and triton (as 



the EDMs of light nuclei in general) are functions of six LECs: 



dsne = 0.88 dn - 0.047 



0.15 1^ + 0.28 |f + 0.01 F^Ci 

IT TV 



0.02 F^C2 efm (59) 



and 



^0 



9i 



rf3H = -OmO dn + 0.90 dp + 0.15|f -0.28|f + 0.01F;Ci -0.02F;C2 efm , (60) 



where ^0 applies for 6, qCEDM, and xl; 9i for qCEDM and xl; dn,p for qEDM and x^'i 
and Ci^2 for xl only. 

Only in the case of the qEDM do we expect the tri-nucleon EDMs to be dominated by 
the nucleon EDMs. Not surprisingly, the helion (triton) EDM should be approximately 
equal to the neutron (proton) EDM. The nucleon EDM for the dimension-six sources was 
calculated in Ref. [26] and it was found that for qEDM the EDMs were dominated by 



the short-range contributions in Eqs. (\21n and (22). In this case 



rf3He(qEDM) 
rf3H(qEDM) 



0.83 rfo - 0.93 til , 
0.85 Jo + 0.95 (Ji . 



(61) 
(62) 



For the 6 term, on the other hand, the helion and triton EDMs depend at LO only 
on ^0- To check this statement we compare the LO contribution with the contribution 
from the nucleon EDMs. If we assume the neutron and proton EDMs to be saturated 
by their long-range part, that is, the chiral log in Eq. (20), Eq. (24) shows that the 



short-range term is comparable to the pion-exchange contribution. To be on the safe side. 
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it seems better not to neglect the LO nucleon EDMs for the 9 term, even though the 
power counting tells us it should be subleading; then 

c/3Hc(^) = 0.83 Jo - 0.93 Ji -0.15 l^efm , (63) 

dmiO) = 0.85 Jo + 0.95 Ji + 0.15 l^efm . (64) 

-TV 

This argument holds equally well for the qCEDM, except that now also gi contributes: 
j3He(qCEDM) = O.83J0-O.93J1- (0.15 1^ + 0.28 If )efm, (65) 

\ F-K Ft, J 

j3H(qCEDM) = 0.85 Jo + 0.95 Ji+ (0.15 If -0.28 If ) efm . (66) 

\ Ft, Ft, J 

Finally in the case of we expect the tri-nucleon EDM to consist of (70,1 pion exchange, 
insertions of the short-range NN interactions, and the contributions from the nucleon 
EDMs. Similarly to the qEDM, the nucleon EDM from is dominated by short-range 
contributions. All six LECs contribute: 



rf3He(xI) = 0.83c;o-0.93Ji- (0.15 |f + 0.28 |f + 0.01 F^Ci- 0.02^^3^2 )efm 



(67) 



ci3H(xI) = 0.85Jo + 0.95Ji+ ( 0.15 |f- 0.28 |f + 0.01 F^Ci- 0.02 F^Csjefm 



(6^ 



From Eq. (|17)) we infer F^Jo,i/^o = 0(eF^/m^) ~ 2 efm and F^Ci.s/^o = 0{F^/ml) ~ 
2. We see again that the ^0,1 coefficients are somewhat smaller than expected; moreover, 
the short-range NN interactions might contribute even less. However, one should keep 
in mind the large uncertainty in the Ci^2 coefficients, and that these dimensional-analysis 
estimates could easily be offset by dimensionless factors in the LECs. 



7 Discussion and Conclusions 

Historically, hadronic EDMs have mostly been discussed in the framework of a one-boson- 
exchange model. It is assumed that P- and T-violation is propagated by pions which are 
parametrized by three fl" non-derivative interactions. In our notation, 

C = -^Nt -TzN - ^Nn^N - ^Nt3TT3N (69) 

jy T, ^ TT ^ IT 

(in the nuclear physics literature, where chiral symmetry and power counting are not 
emphasized, the coefficients are normally defined without Ft,). Hadronic EDMs are calcu- 
lated as a function of these three parameters. In some cases the effects of heavier bosons 
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Table 1: Dependence of the EDMs of the neutron, proton, deuteron, hehon, and triton on 
the six relevant low-energy constants. A "-" denotes that the LEG does not contribute 
in a model-independent way to the EDM at leading order. Values are for the Argonne 
vl8 potential; for the potential- model dependence of the results, see text. 



LEG 


do 


di 


{go/F^)efm 


{gi/F^)efm 


(F^COefm 


(F3C2)efm 


dn 


1 


-1 










dp 


1 


1 










dm 


2 





0.0002 - 0.07/3i 


-0.19 






dsHe 


0.83 


-0.93 


-0.15 


-0.28 


-0.01 


0.02 


dsn 


0.85 


0.95 


0.15 


-0.28 


0.01 


-0.02 



are included as well. In this work we argue that this model is oversimplified. There is 
a priori no reason not to include photon-nucleon and short-range nucleon-nucleon 
interactions at low energies. By studying the chiral properties of the fundamental 
sources of dimension up to six at the QGD scale, it is possible to construct a model- 
independent hadronic Lagrangian with a definite hierarchy between the different 
hadronic interactions. It is found that the one-pion-exchange model with three LEGs is 
not appropriate for any of these sources, and in general there are six hadronic 
interactions that determine the EDMs of light nuclei. Two of those are in the OPE model 
as well — go and gi — and the other four are additional interactions that need to be consid- 
ered when determining hadronic EDMs. The g2 interaction is not relevant at LO for any 
of the fundamental sources. The other four necessary LEGs are the isoscalar and isovector 
components of the neutron and proton EDMs and two isoscalar NN interactions of 
short range. The isovector fl" NN interactions come in at higher order for all sources. 

We therefore propose that nuclear EDMs be analyzed on the basis of these six LEGs. 
In the previous sections we discussed EDMs of light nuclei, providing specific examples in 
the form of the deuteron, helion, and triton. In Table [T] the dependence of these various 
EDMs on the six LEGs is summarized. From the table it is clear that using the OPE 
model gives an oversimplified view. At least six observables are required to identify the 
six LEGs. If other light nuclei become the target of experimental investigation, their 
EDMs can be calculated along similar lines at the cost of larger computer resources. We 
hope that EDMs of heavier systems can be also expressed in terms of these six LEGs. 
However, in these cases there could be significant enhancement factors for the potential 
contribution [50J, making important otherwise subleading terms in the potential [31], such 
as the third non-derivative pion-nucleon coupling ^2- 

Once (a subset of) the LEGs are determined it is possible to learn something about the 
more fundamental sources at the QGD scale. In Table |2] we list for the different 
sources the expected orders of magnitude of the neutron EDM, dn, and ratios between the 
other EDMs considered here and dn- Although some care is needed when using this table 
— as we have discussed, the numbers found earlier are not always exactly of the expected 
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Table 2: Expected orders of magnitude for the neutron EDM (in units of c/Mqcd), the 
ratio of proton-to-neutron EDMs, the ratio of deuteron-to-neutron EDMs, the ratio of 
hehon-to-neutron EDMs, and the ratio of triton-to-neutron EDMs, for the 6 term and the 
three dimension-six sources. Q stands for the low-energy scales F^^, m^r, and 7. 



Source 



9 



qCEDM 



qEDM 




0{1) 


0{l) 


0{1) 


0{l) 


0{1) 


0{l) 


Oil) 


Oil) 



Mqcd dn/e 
dp/ dn 
d^u/ dn 

d-iUe/dn 

dsn/ dn 



09 



^"qcd 



Oil) 
Oil) 



o 



^QCD 



Oil) 



o 
o 
o 



^QCD 



*^QCD 



size — it does allow some qualitative statements, even if less than six measurements are 
available. 

The simplest scenario is the one where is dominated by the qEDM, in which case 



all light nuclear EDMs are essentially given by two LECs only: do and di. (See Eqs. (21) 



(22), (50), (61), and (62).) A measurement of the proton and neutron EDMs would make 
deuteron and tri-nucleon EDMs testable predictions, 

d2ti ~ dn + dp , (70) 
rfsHe + dsn ~ 0.84(d„ + dp) , (71) 
dsHe - rfsH ~ 0.94((i„ - dp) . (72) 

The nucleon Schiff moments and the deuteron magnetic quadrupole moment (MQM) 
depend on other LECs [2Sl EB] and cannot be predicted. For light nuclei the effects of 
the potential from the qEDM are suppressed compared to the nucleon EDMs [3^ . 
although enhancements could make them more relevant for heavier nuclei. 

from 9 and qCEDM manifests itself in EDMs of light nuclei that differ significantly 
from the EDMs of their constituents. For both sources, the EDMs we calculated depend 
at LO on four of the six LECs — Qq, Qi, d^, and di — but in different ways. For the 



qCEDM, the distinguishing feature is that the deuteron EDM (49) is expected to be 
significantly larger than the isoscalar nucleon EDM, thanks to gi. Thus, a measurement 
of nucleon and deuteron EDMs could be sufficient to qualitatively pinpoint, or exclude, 
qCEDM as a dominant source, and to fix the values of (io,i and gi. Then, the isoscalar 
combination of helion and triton EDMs, da^e + da^, which in LO only depends on gi, 
becomes a falsifiable prediction of the theory, 

dsue + dsu ~ 3d2ii . (73) 

If, to be on the safe side, we keep some subleading terms (the nucleon EDMs) as we did 



in Eqs. (65) and (66), then we get (including the subleading term 2dQ in Eq. (49)) an 
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additional — 2.16((i„ + dp) on the right-hand side of Eq. (73). Furthermore, go can then 



be extracted from dsHe — ^sh (see Eqs. (65) and (66)), leading to testable predictions for 
other observables. 

In contrast, for the Standard Model 6 term we do not expect the deuteron EDM to be 
significantly different from twice the isoscalar nucleon EDM. Although the deuteron EDM 



(52) formally depends on the isoscalar nucleon EDM and on the pion-nucleon couplings 



go and gi, the results of Section [5] show that the pion-exchange contribution is likely only 



10% of the nucleon EDM. On the other hand, the EDMs of ^He and ^H, Eqs. (63) 



and (64), are dominated by although they receive important contributions from the 
neutron and proton EDMs. In particular, we expect the isovector combination dan^ ~ da^, 
which is sensitive to to differ from the isovector nucleon EDM di, while the isoscalar 
combination dsHe + da-^ should be close to 2do: 

da-ae + dan ~ 0.84(rf„ + dp) , (74) 
da-He - dan ^ 0.94(rf„ - dp) . (75) 

The experimental observation of these relations in nucleon, deuteron, helion, and triton 
EDM experiments would qualitatively indicate the 6 term as the main source responsible 
for fl". Quantitatively, the measurement of nucleon, helion, and triton EDM allows 
extraction of the coupling go, which then can be used to provide testable predictions of 
other fl" observables, like the proton Schiff moment p6l 120] or the deuteron MQM |28] . 
which are not sensitive to the nucleon EDMs. 

Finally, in the case of the xl sources the analysis is in principle most complicated, due 



to the appearance of all six LECs. Like for 6, the deuteron EDM (51), although formally 



dependent on gi at LO, is probably dominated by do- The tri-nucleon EDMs (67) and 



(68) formally depend on all six LECs, but they are again possibly dominated by do and 
di. It might thus be difficult to separate the sources from qEDM. For less dilute, but 
still light, systems we expect different results. For these systems, in the case of qEDM the 
EDMs are still dominated by (Jo,i, but for sources we expect the contributions from the 
fl" potential to be more significant, implying that measurements on these systems might 
separate xl sources from qEDM. Of course, more extensive calculations are necessary to 
verify this claim. 

In conclusion, we have argued that an experimental program to measure light nuclear 
EDMs could offer valuable information on yet undiscovered sources of parity and time- 
reversal violation. Our case is based on some crucial, but relatively general assumptions, 
such as the validity of the Standard Model with its minimal particle content at the elec- 
troweak scale, and the naturalness of interaction strengths. Elsewhere [28] it has already 
been pointed out — basically on the basis of dimensional analysis — that sensitivity to the 
deuteron EDM at the level hoped for in storage ring experiments [6] would probe scales 
where new physics is expected. A similar analysis holds for our tri-nucleon results, Eqs. 



(59) and (60). But our results here go beyond dimensional analysis and suggest that, at 
least for the lightest nuclei, the contribution of the neutron and proton EDMs are more 
important than expected by simple power counting. For all sources, they compete with, 
when they do not dominate, the effects of the potential. For this reason, other 
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observables insensitive to the nucleon EDMs, for example higher electromagnetic mo- 
ments, could provide important complementary information and a cleaner way to extract 
pion-nucleon and nucleon-nucleon couplings. Additionally, it would be interesting if 
EDMs of heavier systems could be recast in terms of our EFT approach. 
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A Potential in Coordinate Space 



9i9a 
"2F2 

1 



+ r^^)\ (^{i) 



In configuration space, the LO potential of Sect. |3.2| is given by 
%(r) = -^r»-r(^)(<?«-<?(^))-(v.f/(r)) 

) j _ ^0)) + [r^ _ ^0) j (^W + . (^v,t/(r)) 

+ ^ [Ci + C2T» ■ r(^)] (a« - ■ (V.5(^)(f )) , (76) 
where r = Xi — xj is the relative position of the two interacting nucleons and 

U(r) = [2exp (— m^±r) -|- exp (— m^or)] , (77) 

127rr 

which reduces to the usual Yukawa function U{r) = exp(— mjrr)/47rr when, at LO, we 
ignore the pion mass difference. 



28 



Analogously, the NNLO potential of Sect. |4.2| becomes 

go9A 



%(f,V„Vx 



2F2 



3F2 
X (<?« - 



2171-^ 



Vr rU{r)] + {VrW{r 



U{r)Vx (VrVyU{r))V 



where X = {xi + Xj)/2 and 

W{r) 



[exp (— m^rir) — exp {—m^or)] 



which is entirely a consequence of isospin breaking. 



(7J 



(79) 



B Fourier Transform of the Currents 

To evaluate the matrix elements in Section [5] we need to transform the currents to config- 
uration space. We follow Ref. [M] and transform with respect to the nucleon momenta 
but not with respect to the photon momentum. In the most general case 



J (xj, Xj , Xj, Xj,q) 



3 ^-ipi -Xi ^-ipj -Xj ^ip^ -xl ^ip. -x. 



(27r)3 j (27r)3 j (27r)3 j {2Ti f 
{2iTf6'-^\pi+Pj-p- -pj - q)J°{Pi,pi,Pj,Pj,q 



Introducing the relative configuration-space coordinates X 
{xi + Xj)/2, and X' = {x- + Xj)/2, we rewrite this as 

d^R r d^K r d^k 



;o) 



j\f,f',X,X',q: 



(27r)3 J (27r)3 J (27r)3 



^iPtiX-X') -iK-{f~f') 



^-hHr+r')jO(^^^j:^K,Pt) . (81) 



The currents we need (the third currents in Eqs. (36) and (38)) depend on q and k 
only, such that the expression can be simplified to 



jO(f,X,g: 



-iq-X 



d^A; 
(2^ 



e-^^-V°(g,A;) . 



(82) 
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The Fourier transforms can be done and we find for the required currents 

TT 



xe 

F3 



W{q,f), (83) 



xe 



-iq-X 



l^(g,f), (84) 



in terms of the function 



W{q, r' 



8n 



da exp 



i^(l -2a) 



exp[—r{ml + g^a(l — a))^/^] 
[m2 +g2a(l -a)]V2 



^5) 



Before continuing it is convenient to look at the inverse Fourier transform of the current 

= Ze-iq - D + 0{f) , 



^6) 



where Ze is the total charge and D is the EDM operator used in Sections |5j and [6j An 
easy way to extract the EDM operator is by using 



D = ilimVqj\q) . 



g-s>0 



^7) 



As an example we consider the EDM operator coming from J^^{f,X,q). From Eq. 



(87) we read off 

— * 



"^iegAgc 



X 



i?2 

TT 

i 
2 



(a^« + a^0))+X(a«-a(^))-(2V. 



F2 



X 



TTItt \xi Xj I 



where we used limg^o VqVr(g, r) = 0. This is Eq. (43). Following similar steps we obtain 



Eq. disl) from J'^j,^{f,X,q 
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